The boson-fermion model with on-site Coulomb repulsion between fermions 
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The boson-fermion model, describing a mixture of itinerant electrons hybridizing with tightly 
bound electron pairs represented as hard-core bosons, is here generalized with the inclusion of a 
term describing on-site Coulomb repulsion between fermions with opposite spins. Within the general 
framework of the Dynamical Mean-Field Theory, it is shown that around the symmetric limit of 
the model this interaction strongly competes with the local boson-fermion exchange mechanism, 
smoothly driving the system from a pseudogap phase with poor conducting properties to a metallic 
regime characterized by a substantial reduction of the fermionic density. On the other hand, if one 
starts from correlated fermions described in terms of the one-band Hubbard model, the introduction 
in the half-filled insulating phase of a coupling with hard-core bosons leads to the disappearance of 
the correlation gap, with a consequent smooth crossover to a metallic state. 

PACS numbers: 71.10.Fd, 71.30.+h, 71.10.Fd, 74.20.Mn 
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A large variety of phenomena in condensed matter 
physics has been studied in the past in terms of interact- 
ing models with coupled bosonic and fermionic degrees 
of freedom [[tj . Their relevance has proved to be of par- 
ticular evidence in the analysis of the elect ron-phonon 
problem and of the related issue of the polaron and bipo- 
laron formation. In this context, it was proposed some 
years ago that for a system of itinerant electrons in- 
teracting with local lattice deformations, the crossover 
regime between adiabatic and non-adiabatic behavior can 
be described by a model where tightly bound electron 
pairs (hard-core bosons) of polaronic origin coexist with 
quasifree electrons (fermions), with an exchange coupling 
assumed between them by which bosons can decay into 
pairs of itinerant fermions and vice-versa. It has sub- 
sequently been suggested || that this Boson-Fermion 
Model (BFM) can provide a possible scenario for high- 
T c superconductivity, according to the hypothesis that 
the fermionic degrees of freedom describe holes confined 
in the copper-oxygen planes, while the bosonic ones are 
associated with bipolarons which form in the highly po- 
larizable dielectric layers sandwiching the Cu02 planes. 
Since then, the model has been deeply investigated, in 
particular by Ranninger and coworkers [@-f| , with a spe- 
cial attention to the description of the pseudogap phase 
characterizing the normal state of the underdoped high- 
T c copper oxides. 

In this paper we consider an extension of the BFM 
in which fermions, besides exchanging with bosons, are 
assumed to also experience an on-site Coulomb repul- 
sion. This interaction, tending to prevent two fermions 
with opposite spin from occupying the same site, strongly 
competes with the boson-fermion exchange mechanism, 
giving rise to interesting effects in the spectral as well as 
in the transport properties. The model, whose fermionic 
part goes back in the one-band Hubbard model in the 
limit of vanishing boson-fermion coupling, is solved here 



within the framework of the Dynamical Mean-Field The- 
ory (DMFT) by applying to the related self- 
consistent impurity problem the so-called Non-Crossing 
Approximation (NCA) @. We recall that DMFT-based 
approaches have been intensively used in the last years in 
the related context of clcctron-phonon systems described 
by the Holstein model Jl3|,[l4| or extensions of it |D| , as 
well as in the case of boson-fermion coupled systems with 
boson dispersion |l6| . By means of various analytical and 
numerical techniques, a special attention has also been 
devoted to the role played by the on-site Coulomb repul- 
sion in the electron-phonon problem and, in particular, 
on the polaron and bipolaron formation. This has been 
done essentially referring to a generalized version of the 
Holstein model which includes an Hubbard interaction 
term (the so called Holstein-Hubbard model) (l7) . 

The hamiltonian for the BFM with Coulomb repulsion 
between fermions is given by 

H = eoY^ cL c i* ~ 1 4a c j<r + U J2 4^4x^1 

+ EoY, b l b i + gY,^iCiT + 4Ai b ^ ■ W 



Here c^J denote annihilation (creation) operators for 
electrons with spin a at site i and denote hard-core 
bosonic operators describing tightly bound localized elec- 
tron pairs. The site energies for fermions and bosons are 
expressed as £q = D — li and Eq = — 2/i, respec- 
tively, where the chemical potential li is assumed to be 
common to the two kinds of carriers (up to a factor 2 
for the bosons) in order to guarantee charge conserva- 
tion. Finally, g denotes the boson-fermion pair-exchange 
coupling constant, while t, D, and U are respectively 
the bare hopping integral, the bare half bandwidth and 
the on-site Coulomb repulsion, all referred to fermions. 
We see that in the limit of vanishing Coulomb repul- 



sion the model goes back in the boson-fermion model for 
high-T c superconductors analyzed in Refs. On the 

other hand, in the limit of zero boson-fermion coupling 
the fermionic part of the model defined by the Hamilto- 
nian (^) coincides with the standard one-band Hubbard 
model. 

It is worth pointing out that the assumption of an 
on-site Coulomb interaction affecting only the fermionic 
degrees of freedom is motivated by the existence of 
anisotropic systems, such as the previously mentioned 
high-T c copper oxides, where there is experimental evi- 
dence that bosons and fermions reside in separated lattice 
regions where Coulomb effects are in general significantly 
different. In this respect, the site indices in the phe- 
nomenological Hamiltonian (1) would better be thought 
as representative of molecular units composed, for in- 
stance, by an atom on which the correlated fermions sit 
and two adjacent ones forming a diatomic molecule hous- 
ing the bosonic tightly bound electron pairs. Given this 
picture, no direct interference occurs between the on-site 
fermionic Coulomb repulsion and the bosonic formation, 
the only effect being a shift at higher energies of the 
fermionic level with respect to the bosonic one, leading 
to a relative change in occupation of the two species. 
The assumption of a lattice seen as a collection of effec- 
tive sites is of course made to simplify the analytical as 
well as the numerical calculations, and a more accurate 
(though technically more complicate) way of analyzing 
the role of the Coulomb repulsion would require the so- 
lution of the model on a true bipartite lattice. 

Given the local character of the interaction terms in 
the Hamiltonian (Q), important insights on the physics 
of the model can be obtained from its atomic limit. The 
eigenstates and eigenvalues of H are in this case given by 
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where the dot indicates the presence of a boson and we 
have defined 
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• = 1. When C/ = and the total particle 
density n = rip + 2ns is equal (or close) to 2, the physics 
underlying the model is fully determined by the charge 



exchange interaction. As the temperature is decreased, 
the spectral weight of the pole in the Green's function 
associated with the non-bonding single-particle state (at 
energy £0) is gradually transferred to the poles (at ener- 
gies E4 — £0 and E5 — £0) associated with the bonding 
and anti-bonding two-particle states |4 > and |5 > p8[ . 
The corresponding effect seen when itinerancy is taken 
into account, is a depletion of the density of states in the 
region around the chemical potential. This gives rise to 
the opening a pseudogap which, according to the mech- 
anism explained above, is thus induced by local pairing 
among the electrons. 

When the effect of the on-site Coulomb repulsion 
among fermions is introduced, the structure of the atomic 
Green's function becomes more complex, with the num- 
ber of poles increasing from 3 to 6. The transfer of spec- 
tral weight from single-particle to two-particle resonances 
as T is decreased, becomes less and less pronounced as 
higher values of U are considered. As expected, this is 
a consequence of the fact that a finite U tends to forbid 
double fcrmion occupations of the same site, thus inhibit- 
ing the boson-fermion charge exchange mechanism. More 
precisely, what one can see from the atomic limit solution 
is that, upon increasing U , u 2 tends to zero, v 2 tends to 
1, and thus the two states |4 > and |5 > tend to become 
pure one-boson and two-fermion states, respectively, with 
a separation in energy of the order of U. As we will see 
in the following, for the full model around the symmet- 
ric limit this has the simultaneous effect of hindering the 
opening of the pseudogap and making the bosonic site 
occupation considerably higher than the fermionic one. 

The effect of the electron itinerancy will now be taken 
into account within the general framework of the Dynam- 
ical Mean-Field Theory (DMFT) In the last years 
this approach has proved to be of fundamental impor- 
tance in condensed matter physics and has been success- 
fully applied to a variety of many-body models. Within 
DMFT an interacting system is seen as a purely local 
system coupled to a Weiss "mean field" to be determined 
self-consistently, produced by the neighboring sites. This 
approach leads to the freezing of all spatial fluctuations, 
just like the usual classical mean-field theory, but has 
the advantage of fully retaining the quantum temporal 
fluctuations of the original problem (hence the adjective 
"dynamical" ) . Indeed, the mean field is here a function 
of time, instead of being a pure number, associated with 
the probability amplitude of creating a particle on the 
impurity center at a given time and destroying back in 
the external bath at a later time, and vice- versa. Like any 
mean field approximation, this approach becomes more 
and more accurate as the lattice coordination number in- 
creases, becoming exact in the limit of infinite dimension- 
ality. This basic structure makes the application of the 
DMFT particularly suited to models, such as that defined 
by the Hamiltonian (Q), where the interaction terms are 
local. A method for the inclusion of non-local corrections 
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to DMFT, known as Dynamical Cluster Approximation, 
has recently been developed pl| and applied to the two- 
dimensional Hubbard model EG] , by mapping the lattice 
problem onto a periodic cluster embedded in an external 
host. 

Since we consider here a model where correlated elec- 
trons are coupled to dispersionless bosons, in the DMFT 
self-consistent procedure the bosonic channel for the 
Green's function remains frozen, in the sense that the 
related bare Weiss self-energy at each iteration is un- 
renormalized and fixed at its initial non-interacting ex- 
pression. This point, discussed at length in Ref. |l(|, im- 
plies that the general dynamical mean-field framework in 
which the model defined by the Hamiltonian ([!]) is inves- 
tigated, is essentially the same as that developed for the 
Hubbard model in Refs. |7j-|9||, the effect of the boson- 
fermion exchange coupling manifesting itself only in the 
solution of the impurity problem. In the following we 
give a brief summary of the DMFT technique, along the 
general lines of the diagramatic approach presented in 
Refs. ||||. The reader is referred to these papers, as well 
as to Refs. [|],[l(|[ll]] , for all technical details not reported 
here for brevity. 

The approach starts from the demonstration [ pl|j22| 
that in the limit of infinite dimensionality only the lo- 
cal dynamics remains nontrivial, in the sense that the 
only non-vanishing contributions to the self-energy are 
the site-diagonal ones. As a consequence, the Dyson 
equation in real space takes the form 



Gij(iuj n ) = G° tj (iuj n ) + y] Gji (w n )S(w w )Gi!j (iuj n ) 



(4) 



When this equation is referred to the local propagator 
Gii, one can first consider the self-energy contributions 
coming from all sites j other than the site i. The resum- 
mation of these contributions lead to the definition of a 
modified propagator, that we denote by Gw (the Weiss 
field), that is then used in the full Dyson equation for 
Gu, reintroducing the missing self-energy contributions 
that involve the same site i ||. This procedure leads to 
the following Dyson equation for Gu = Gi mp 



G imp {iuj n ) = [G w (iu n ) - E(mj„)] 



(5) 



which proves the equivalence of the lattice problem to a 
model of an Anderson impurity sclf-consistently embed- 
ded in a medium specified by the Weiss Green's func- 
tion Gw- Given for the latter the structure Gw(iu n ) = 
[iu> n — Eq — Eiy(jw n )] _1 , the above equation implies that 
the impurity Green's function can be expressed in the 
form 



Gimp{i^n) 



1 



iUJ n - EQ - XV(i^n) - ^int(i0Jn) 



, (6) 



where we have put Ej„t = E. We thus see that the to- 
tal self-energy is written as the sum of two momentum- 
independent contributions Ej„t and E^y, associated re- 
spectively with the on-site interactions (the boson- 
fermion exchange and the Coulomb repulsion among 
fcrmions) and with the hybridization of the impurity cen- 
ter with the medium. 

On the other hand, due to the locality of the self- 
energy, the lattice Green's function in k-space is given 
by 



G; Qt (£ k , iu n ) 



1 



£a — £k — Ej nt (io; ra ) 



(7) 



where e k denotes the bare electron dispersion. The 
self-consistency condition which iteratively leads to the 
renormalization of the Weiss self-energy, thus closing the 
DMFT algorythm, is then simply obtained by equating 
the impurity Green's function (^) to the integral over k- 
space of the lattice Green's function, i.e. Gi mp [iuj n ) — 
N~ x Giat(£\t, ioJ n )- Upon replacement of the integra- 
tion over k by an integration over energy, this condition 
takes the form 



G imp (iu> n ) = I de 



P(e) 



£o 



E, n t (iuj n ) 



(8) 



where p is the free density of states associated with the 
particular lattice chosen. 

The whole DMFT procedure can thus be summarized 
as follows. One starts with a guess for the Weiss Green's 
function Gw and solve the corresponding impurity prob- 
lem, obtaining Gi rnp (here this is done by means of the so- 
called Non-Crossing Approximation). The resulting local 
self-energy Ej nt = E, evaluated from Eq.(||), is then used 
to find a new impurity Green's function by imposing the 
self-consistency condition (||). By subtracting off from 
the latter the local self-energy, again by use of Eq.(||), 
one gets a new Weiss field which is used to reinitialize 
the process, which is then iterated until convergence is 
reached. 

Let us also notice that the explicit form taken by the 
self-consistency condition depends on the particular lat- 
tice on which the DMFT equations are solved. In the 
following wc assume to refer to a Bethe lattice with co- 
ordination number z —> oo and nearest neighbor hopping 
tij = t/^/z, for which the electronic bare density of states 
is the semi-circular function p(e) — (l/27rf 2 )\/4i 2 — e 2 
(|e| < It) of width 2D = At. This assumption, widely 



used in the context of the DMFT |j-[jj|^,|2J allows to 
solve analytically eq. (jsl), leading to 



T,w(iw n ) = t 2 Gim P (iu) n ) 



(9) 



or, in terms of the corresponding imaginary parts, 

A(w) = i 2 A f (uj) (10) 



(we use the definition A(uj) — — 2ImSvi/(w + iff), with a 
similar equation relating Ap and Gi mp ). 

The single-site effective dynamics described above is 
naturally investigated by means of functional methods 
based on the use of an imaginary-time effective action for 
the fcrmionic degrees of freedom defined at the impurity 
center |7|Jl0|. An equivalent picture can be introduced 
by considering a related Hamiltonian formulation where 
the Weiss field is described in terms of auxiliary degrees 
of freedom. In the case of the BFM, we have an effective 
Hamiltonian of the form 

a 

+ Eotfb + g [ c\c\b + 6 + c;c T ] 

+ ^2w k dl a d ka + ^2 v k [d\ a c a + c\d ka } , (11) 

k , a k,{7 

where c& and denote the original fermion and bo- 
son operators, respectively, at the selected site and d k ^ 
are the auxiliary fermionic operators associated with the 
Weiss field, taking into account the dynamics occurring 
at all the other sites. Their energy spectrum w k and their 
coupling v k to the impurity electrons are of course not 
known a priori but must be determined self-consistcntly. 
Identifying the Weiss Green's function with the propaga- 
tor for the auxiliary variables, one gets 



Gw(iu> n ) = iuj n - g - 



lUJ n - W k 



(12) 



with the sum at the right-hand side representing the self- 
energy Yjw . We thus see that the Hamiltonian descrip- 
tion is fully equivalent to the functional one provided 
that the parameters v k and w k are related by the fol- 
lowing equation to the imaginary part A(cl>) of the Weiss 
self-energy appearing in the effective action: 



A(cj) = 2ir v k S(u> — w k ) 



(13) 



Depending on the technique adopted to solve the impu- 
rity problem, the Hamiltonian formulation can prove to 
be useful in practical calculations, besided providing a 
description probably closer to our physical intuition. 

The self-consistent single-impurity Anderson problem 
defined by the effective Hamiltonian ( |TT| ) is solved here 
within the so-called Non Crossing Approximation (NCA) 
fL2f , along the lines presented and discussed in Ref. [^). 
Within the DMFT framework, NCA has been applied in 
the last years to several other models, such as the one- 
band Hubbard model |2q| , a multiband Hubbard model 
for perovskites |2(J , the Anderson lattice model with cor- 
related conduction electrons J2j|, and the Kondo lattice 
model with correlated conduction electrons [Q. In a 
very recent paper this technique has also been employed 



to analyze the electronic structure of the ruthenate alloy 
series Ca2- a; Sr ;r Ru04 p7| . 

All the results presented in the following refer to the 
case of a total particle density n — rip + 2ns equal to 
2, with np and Hb being respectively computed from 
the fermionic and the bosonic impurity propagators. All 
energies are measured in units of the bare bandwidth 
2D. In order to emphasize the effects of the competi- 
tion between the boson-fermion coupling and the on-site 
Coulomb repulsion, we will mainly consider parameter 
regimes where g and U are comparable. We want also 
to point out that in these regimes the application of the 
conventional NCA approach is expected to give reliable 
results. This has been checked resorting to a validity cri- 
terion which is discussed at length in ref. || and which, 
for brevity, will not be repeated here. Finally, we have 
also verified that in the limit of vanishing boson-fermion 
coupling, the results for the Hubbard model obtained 
in Ref. |§5|, where the same DMFT+NCA algorythm is 
used, are correctly reproduced. 
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FIG. 1. Fermionic density of states for T = 0.06, g = 0.4 
and several values of U. In the case U = the site energies 
satisfy the symmetric limit condition Eo = £o = 0. 

We start by discussing how the inclusion of U affects 
the pseudogap phase as described by the BFM studied in 
Refs. In Fig. |] we have plotted the fermionic den- 

sity of states for several values of the on-site repulsion 
(the energies are measured with respect to the chemi- 
cal potential). We start at U — from the fully sym- 



metric case Eq = Eq = (top panel), which implies 
Uf = 2riB = 1, and then consider finite values of U while 
keeping the bosonic bare energy Ab fixed (the finite-C/ 
symmetric limit conditions Eq = 2eq + U = are con- 
sequently not satisfied). We see that the widely open 
pseudogap showing up for U = is gradually filled as 
increasing values of U are considered. At the same time 
an Hubbard correlation gap of the order of U develops 
well beyond the chemical potential, thus without sub- 
stantially affecting the low-energy physics. As already 
pointed out in the discussion of the atomic limit, this 
is due to the fact that the on-site repulsion, tending to 
inhibit the boson-fermion exchange mechanism, pushes 
the system in a regime where the coupling between the 
two components is reduced and, for the parameter choice 
considered here, the almost pure bosonic level remains 
pinned below the fermionic ones (these latter being split- 
ted by U). Starting from the U = fully symmetric 
case considered in the upper panel, this also implies that 
much of the fermionic spectral weight is transferred at 
higher and higher energies as U is increased, this lead- 
ing to a substantial reduction of the fermionic occupation 
with respect to the bosonic one. We thus see that for the 
system is energetically more favourable to move away 
from the "insulating" half-filled case, creating holes in 
the fermionic configuration and thus allowing a smooth 
transition towards a metallic behavior. 
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FIG. 2. Real part of the optical conductivity for the same 
parameters as in Fig^ and several values of U. 

These considerations find further support in the be- 
havior of the real part of the optical conductivity, which 
in the limit of infinite dimensions becomes directly con- 
nected to the fermionic single-particle spectral function 
Ap(ek,uj) according to the relation 

a(oj) — 7r de p(e) I dcu' A F (e, ui')A F (e, lj + w') x 



(14) 



x — \np(uj') — npiuj + u/)l 



ported in Fig.0 we see that a situation typical of that 
expected for an insulator is found for U — ► 0, that is, 
we obtain a small d.c. conductivity er(0) with a peak at 
uj ~ 2g reflecting transitions between the bonding and 
the antibonding resonances exhibited by the density of 
states (see Fig.|l|). This peak gradually disappears as U 
is increased, with a concomitant Drudc-like accumulation 
of weight al low uj, clearly associated with the establish- 
ment of metallic properties. 
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(here the sum over momenta has again been expressed 
as an energy integration over p(e)). From the results re- 



FIG. 3. Fermionic density of states in the symmetric limit 
(So = 2e + U = 0) for T = 0.04, U = 1 and several values of 
9- 

We have seen so far that if one starts from the symmet- 
ric limit of the BFM, characterized by a charge- exchange 
induced pseudogap in the excitation spectrum, then the 
introduction of the on-site Coulomb repulsion between 
fcrmions leads to a transition from an insulating state to 
a metallic state with lower carrier density. Recalling that 
in the limit .9 = the fermionic part of the model defined 
by the Hamiltonian (Q) becomes equal to the Hubbard 
model, we now reverse the point of view, analyzing how 
the physics of the Hubbard model in the half-filled insu- 
lating phase is modified by the inclusion of the boson- 
fermion coupling. To this purpose, we assume to start 
from a parameter regime where, in the absence of g, the 
correlation gap is fully openin a half-filled configuration. 
Then we turn on the effect of g under the assumption 



that the bosonic energy level is superimposed just at the 
middle of the unperturbed fermionic band, in such a way 
to ensure a particle density equally distributed among 
fermions and bosons (hf = 2ns = I). As previously 
mentioned, this is realized when the symmetric limit con- 
ditions Eq — and 2so + U = are simultaneously satis- 
fied. The corresponding behavior of the fermionic density 
of states and of the real part of the optical conductivity 
is illustrated in FigsJ| and ||, respectively. 
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FIG. 4. Real part of the optical conductivity for the same 
parameters as in FigJ^. 

We see that, as g is increased from zero, the Hubbard 
gap is gradually closed as a consequence of the fact that 
the boson-fermion coupling tends to favor configurations 
with double fermionic occupations. As a result, this 
mechanism tends to create empty sites in the half-filled 
configuration, thus promoting the fermionic itinerancy. 
For suitably high values of g a narrow pseudogap at the 
chemical potential tends to form again, driving back the 
system in a state where the metallic behavior is gradu- 
ally weakened as T is decreased. The substantial change 
in the transport properties induced by the boson-fermion 
coupling is also evident from the temperature behavior 
of the dc resistivity, plotted in FigJ| (due to numerical 
difficulties our NCA algorythm is currently not able to 
provide stable enough results at temperatures lower than 
those considered in the figure). 
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FIG. 5. Temperature dependence of the dc resistivity in 
the symmetric limit for (7 = 1 and several values of g. 
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FIG. 6. Imaginary part of the fermionic self-energy for the 
same parameters as in Fig.tA 



We can see that, while at g = and for low values of 
g the system exhibits the expected semiconductor-like 
behavior, with <t(0) _1 going up as T is decreased, on 
the other hand higher values of g give rise to a metallic- 
like behavior characterized by a decrease of cr(0) _1 as 
T — > 0. Nonetheless, an important qualitative difference 
should be noted between the cases of intermediate and 
high values of g. The development of a dip at the chem- 
ical potential above approximately g = 0.4 pushes the 
system toward a state which, though still showing metal- 
lic properties, is characterized by strong deviations from 
the Fermi liquid behavior. This is particularly evident 



when plotting the imaginary part of the fermionic self- 
energy for the same parameter choices as in Fig|| We 
can see from Fig.^Jthat -ImEj n i(w) exhibits for g = 0.3 
a parabolic minimum at /i (though with lm£(0) m ( 7^ 0) 
which is the manifestation of a regime close to that of a 
normal Fermi liquid. When higher values of g are con- 
sidered, the behavior becomes significantly different, with 
the minimum being replaced by a narrow peak associated 
with a resonant scattering of the fermionic quasiparticles 
with the bosons. 

In conclusion, we have analyzed a two-component 
model where correlated fermions and localized hard-core 
bosons interact through a mechanism by which bosons 
decay in two fermions with opposite spins and vice- versa. 
A special attention has been devoted to how the inter- 
play between the boson-fermion coupling and the on-site 
Coulomb repulsion between fermions affects the spectral 
and the transport properties. Interesting crossovers from 
a pseudogap phase to a metallic regime, and vice-versa, 
have been found, as driven by the competition between 
the above mentioned local interactions. We emphasize 
that the way we have taken into account the effect of the 
on-site Coulomb interaction among fermions is not real- 
istic if the model is applied to the description of high-T c 
superconductors. A more proper way would rather con- 
sist in defining the model on a lattice with two kinds 
of sites, one associated to the dielectric layers and the 
other one to the copper-oxygen planes. The Coulomb 
repulsion should then be introduced to describe correla- 
tion effects among holes in the latter, with the charge 
exchange mechanism now involving two strongly corre- 
lated fermions sitting on neighboring sites, rather than 
on the same site. A study along these lines is planned 
for the near future. 
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